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ABSTRACT 

Analyzing the halo and galaxy catalogs from the Millennium simulations at 
redshifts z = 0, 0.5, 1, we determine the alignment profiles of cluster galaxies 
by measuring the average alignments between the major axes of the pseudo in- 
ertia tensors from all satellites within cluster's virial radius and from only those 
satellites within some smaller radius as a function of the top-hat scale difference. 
The alignment profiles quantify how well the satellite galaxies retain the memory 
of the external tidal fields after merging into their host clusters and how fast 
they lose the initial alignment tendency as the cluster's relaxation proceeds. It is 
found that the alignment profile drops faster at higher redshifts and on smaller 
mass scales. This result is consistent with the picture that the faster merging 
of the satellites and earlier onset of the nonlinear effect inside clusters tend to 
break the preferential alignments of the satellites with the external tidal fields. 
Modeling the alignment profile of cluster galaxies as a power-law of the density 
correlation coefficient that is independent of the power spectrum normalization 
(erg) and demonstrating that the density correlation coefficient varies sensitively 
with the density parameter (Q m ) and neutrino mass fraction (f u ), we suggest 
that the alignment profile of cluster galaxies might be useful for breaking the 
fi m -cr 8 and f u -a 8 degeneracies. 

Subject headings: cosmology:theory — large scale structure of Universe 



1. INTRODUCTION 



The most salient feature of the large scale structure in the Universe is the web-like 
filamentary pattern of the spatial distribution of the galaxies, which is often referred to 
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as the cosmic web (jBond et al.l Il996l ) . Ever since the work of lHawley fc Peebles! (119751 ) 
who proposed testing the statistical significance of the anisotropy in the orientations of 
the galaxies embedded in the large-scale structure, plenty of observational and numerical re- 
searches have been conducted to understand the origin, property and dynamics of the cosmic 
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It is now understood that the anisotropic spatial correlat ions of the gravitational tidal 



fields mould the large scale structure into the cosmic web (e.g.. lAltay et al.ll2006t IWang et al. 



201ll ). which interconnects the Universe on different scal es and transla t es the large-scale 
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Many observational phenomena can attribute to the interconnection between the large- 
scale anisotropy and the small-scale intrinsic alignments. For instance, iGodlowski fc Flin 



( 120101 ) reported the detection of strong correlations between the orientati ons of the galaxy 



groups and the neighbor group distributions on the scale up to 20 /i _1 Mpc. iPaz et al.l (]201lh 
found a strong alignment signal between the projected major axes of the group shapes and 
the surrounding galaxy distribution up to scales of 30 /i _1 Mpc. They noted that this observed 
anisotropy signal becomes larger as the galaxy group mass increases, which is in excellent 
agreement with the prediction of the ACDM cosmology. Their result is also consistent with 
the previous findings that the degree of the alignment of the cluster galax i es depends on the 



velocity dispersion and number of the clu st er galaxies (IPlionis et al.ll2003t lAryal et al.ll2007 



Godlowski et al.l I201CH : iGodlowskil l2011bf ) . ISmargon et al.l ( 1201 ll ) detected a weak signal of 



the alignments between the projected major axes of the neighbor clusters on the scale of 
20ftr 1 Mpc as well strong signal of the alignments between the cluster maior axis and 
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Varela et al.l ( 120121 ) found that the spin axes o f the spiral galaxie s in the vicinity of the voids 
tend to lie on the surfaces of the voids, while IJones et al.l (j2010[ ) showed that the spin axes 
of the spiral galaxies located in the cosmic filaments tend to be aligned with the elongated 
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According to the previous works, the nonlinear evolution of the gravitational tidal fields 
has dual effects. On the one hand it sharpens the cosmic web and enhances the degree of the 
large-scale anisotropy by inducing filamentary merging and accretion. On the other hand 
it diminishes the interconnections between small and large scales by developing secondary 
filaments and local effects. To explain the formation and evolution of the large-scale structure 
in the cosmic web, it is first necessary to understand how the correlation of the tidal fields 
evolves in a scale dependent way. The galaxy clusters are the best target for investigating 
the evolution of the tidal connections between large and small scales. Often forming at the 
intersections of the filaments, the majority of the clusters a re not isolated but i n constant 
interactions with the surrounding large scale structure (e.g., 



Springel et all 120051 ). 



The anisotropic merging of satellites from the surroundings along the major filaments 
leads to the preferential locations of the satellite galaxies near the minor principal axes of 
the external tidal fields (i.e., the directions of the least compression). As the clusters become 
relaxed, however, the internal tidal fields develop into dominance, gradually rearranging the 
satellite distributions and deviating them from the directions of the external tidal fields. 
Those satellites which merged earlier, being located in the inner regions of the clusters, tend 
to become less aligned with the directions of the external tidal fields. In this picture, it 
may be expected that the profiles of the alignments between the spatial distributions of the 
satellites in the inner and the outermost regions of the clusters reflec t how dominant the 
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Here, we determine the alignment profile of cluster galaxies using the numerical data 
from a large N-body simulation and construct a theoretical model for it. The organization of 
this paper is as follows. In §2, the numerical datasets from N-body simulations are described 
and the alignment profiles of the cluster galaxies for the three and two dimensional cases 
are determined. In §3, a theoretical model for the alignment profile of the cluster galaxies 
is presented and fitted to the numerical results. In §4, we discuss the possibility of using 
the alignment profile of cluster galaxies as a cosmological probe. In §5, the results are 
summarized and a final conclusion is reached. 
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2. DATA AND ANALYSIS 

We use the semi-analytic galaxy catalogs as well as the halo catalog^ at z — 0, 0.5 and 1 
from the Millennium Simulation for a ACDM concordanc e cosmology with Q m = 0.25, = 
0.75, h = 0.73, n s = 1.0, <j$ = 0.9 (jSpringel et al.ll2005l ). The dark matter halos and their 
substructures in the catalogs were found with the hel p of the standard friends-of-friends 
(FOF) and SUBFIND algorithms (jSpringel et al.ll200ll ). respectively, while the galaxies in 
the semianalytic catalogs were modeled accor ding to the merger t rees derived from the semi- 
analytic simulations of the galaxy formation (ICroton et al.ll2006l ). 

For the cluster halos with mass larger than 10 14 h~ l M Q at each redshift, we locate 
their satellite galaxies from the semi-analytic galaxy catalog and measure the pseudo inertia 
tensor, (%), defined in terms of the relative positions of the satellite galaxies located within 
cluster's virial radius r vir as 



lij yvir) 



T 



where x r 



x a> i) is the position vector of the a-th satellite relative to the halo center 
XqjI is its magnitude. Diagonalizing (Iij), we find the three eigenvalues and 



and x a 

corresponding unit eigenvectors of the pseudo inertia tensor for each cluster halo. We select 

1 /2 

only those clusters which satisfy {q$/ Qi) < 0.9 where A3 and Ai are the smallest and the 
largest eigenvalues of L 



.), expecting their pseudo inertia tensors to have three distinct 
eigenvectors. For each selected cluster halo, the major axis at the virial radius, e(r vir ), is 
determined as the unit eigenvector corresponding to Q\. 

Using only those satellites located within some radial distance smaller than the cluster's 
virial radius, r < r V i r , we recalculate the pseudo inertia tensor of each cluster's interior region 
measured at r, 



'ij(r) 



£ 

x a <r 



Xct,iXa,j- 



(2) 



Diagonalizing Iij(r) and determining its unit eigenvector corresponding to the largest eigen- 
value, we find the major axis, e(r), at the radial distance r, of the interior region of each 
cluster. The square of the cosine of the angle between the two major axes is now calculated 
as cos 2 ^ = |e(r„j r ) • e(r)| 2 for each selected cluster halo. Figure [1] illustrates the projected 
positions of the satellite galaxies in a selected cluster, showing how the alignment angle, 9, 
between the distributions of the satellites located within two different radial distances, r vir 
and r is measured. Varying the value of r from 0.9r„j r down to some cutoff radial distance, 



1 They are available at http://www.rnpa-garching.mpg.de/millennium 
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r c , defined as the radial distance within which less than five satellites are located, we repeat 
the whole calculation. 

Since we are going to relate the alignment profiles of cluster galaxies to the auto- 
correlations of the tidal fields smoothed on two different top-hat filtering scales (see §3), 
we would like to express cos 2 9 as a function of the top-hat scale difference rather than the 
radial difference r v { r — r. For a given mass M(< r) enclosed by the radial distance r, the 
corresponding top-hat radial scale, R, is obtained as R(< r) = [3M(< r)/(47rp)] 1//3 where 
p is the mean mass density of the Universe. Here, the mass, M{< r ), can be calculated 
th rough integra t ing th e NFW density profile (INavarro et al.lll996l . 119971 ) from to r, as done 



in 



Klvpin et all (119991 ): 



M(< r) 



M„ 
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where r = r/r s and f(r 
M vir = [(4vr/3)p c 200 r 3 



= ln(l + f) — f/(l + f) Here the virial mass M vir is defined as 
where p c is the critical density at z and r s is the scale radius, 
and c is the concen tration parameter defined as c = r vir /r s , depending on M vir and z as 
klvpin et alill999l ) 



124 



1 + z \ IHMb 



-0.084 
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From here on we express cos 2 9 as a function of AR = Rq — R, where Ro and R are 
the top-hat filtering scales corresponding to M vir at z = and M(< r), respectively, which 
can be calculated from r vir and r, by Equations (J3])- ((^J) . The alignment profile of cluster 
galaxies is now defined as (cos 2 9)(AR) — 1/3, where the constant factor 1/3 is subtracted 
from the ensemble average since the value of the ensemble average will be 1/3 when there is 
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To investigate how the alignment profile of cluster galaxies depends on the mass scale, 
we divide the selected clusters at z = into three samples (I, II and III) according to 
their masses and separately measure the alignment profiles from each sample. Table [1] lists 
the redshift, logarithmic mass range and number of the selected clusters belonging to each 
sample in its second, third and fourth columns, respectively. The left panel of Figure |2] plots 
the alignment profiles of cluster galaxies from the subsamples I, II, and III as open triangles, 
diamonds, and squares, respectively. The uncertainties associated with the numerical data 
points represent the Jackknife errors. We divide the simulation volume into eight subvolumes 
(the Jackknife resamples) each of which contains the same number of halos and then measure 
the alignment profile using the halos from each Jackknife resample, separately. The Jackknife 
errors are finally obtained as the one standard deviation scatter among the eight resamples 
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( IWall fc Jenkinal2008[ ). It is worth mentioning here that the alignment profile is plotted as 
a function of AR/R rather than AR (where Rq is the top-hat filtering scale corresponding 
to M vir at z — 0) to make a fair comparison of the results from the three samples which 
must have different values of Rq. As can be seen, there are clear signals of alignments in the 
small AR section, but the alignment profiles drop gradually as AR increases for all three 
cases. The overall amplitude of the alignment profile shows a weak tendency to increase as 
the cluster's mass increases. 

We also measure the alignment profiles of cluster galaxies at z = 0.5 and 1 to investigate 
how the profiles change with redshift. The left panel of Figure [3] compares the three alignment 
profiles at z — 0, 0.5 and 1 (open squares, triangles and diamonds, respectively), where the 
cluster's masses lie in the narrow range of 10 14 < log(M/[/i -1 M ]) < 10 14 ' 2 . As can be seen, 
there is a strong tendency that the alignment profiles at higher redshifts drop faster with 
AR, which implies that at higher redshifts the satellites in the inner regions of the clusters 
tend to be less aligned with the external tidal fields. 

Recalling that what we can readily observe and measure is the two dimensional spatial 
distributions of the cluster galaxies rather than the three dimensional ones, we also measure 
the alignment profiles of cluster galaxies in the two-dimensional projected space. Adopting 
the flat-sky approximation and taking the z-axis as the line-of-sight direction, we project 
the positions of all satellites onto the x-y plane and follow the exactly same procedures 
to determine the two dimensional alignment profiles as (cos 2 8)(R — Rq) — 1/2 where the 
constant term 1/2 represents the case of no alignment for the two dimensional case. The left 
panels of Figures 0] and [5] plot the same as the left panels of Figures [2] and [3j respectively, 
but for the two dimensional case. As can be seen, the overall shapes and behaviors of the 
two dimensional alignment profiles are quite similar to the three dimensional ones except for 
the decrement of the amplitudes. 

It is worth mentioning here that we evaluate the pseudo inertia tensors for the deriva- 
tion of the cluster g alaxy alignment profi le (Eqs.[T]-[2]) rather than the conventional mass- 



weighted ones (e.g., lAllgood et al.l 120061 ). given that the pseudo inertia tensors have been 



found to quantify w ell the anisotropy in the galaxy spatial distribution caused by the lo- 



cal tidal fields (e.g., IShandarin et al.l l2006t iPark fc Led 120071 ). Since the calculation of the 



pseudo inertia tensors does not require any information on the satellite masses, the two di- 
mensional alignment profile of the cluster galaxies is indeed a readily measurable quantity. 
In the following section, we provide physical explanations for the numerical results obtained 
here and present an analytic formula for the alignment profiles of cluster galaxies. 
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A THEORETICAL MODEL 



Assuming that the alignment profiles of satellite galaxies in clusters are directly related 
to the spatial correlations of the tidal fields, we first measure the alignments between the 
minor principal axes of the tidal fields smoothed on two different scales, usi ng the numerical 
data o f the Millennium density fields constructed on 256 3 grids. As done in iLee fc Erdogdu 
( 120071 ). we perform the Fourier transform of the Millennium density field S fx) on each grid 
with the help of the FFT (Fast Fourier Transformation) method ( jPress et al.l 119921 ) and 
calculate the Fourier transform of the tidal tensor smoothed on a given scale Ro as Ty(k) = 
kikjS(h)W(kR )/k 2 where k = (ki) is the wave vector and W(kR Q ) is the top-hat window 
with the filtering radius Rq. The inverse Fourier transform of Tjj(k) gives the real space 
tidal field lij(x) smoothed on the scale of Rq at each pixel point. The unit eigenvector of 
Tjj(x) corresponding to the smallest eigenvalue represents the minor principal axis, u(Rq). 

Repeating the whole process but smoothing the tidal field on some smaller scale R < R , 
we also determine the minor principal axis, u(R), of the tidal field on scale R. The alignment 
between the minor principal axes of the tidal fields smoothed on two different scales, R and 
Ro, is calculated as |u(i?) -u(_R )| 2 at each pixel. Varying the value of R, we repeat the whole 
calculation, average the alignments over 256 3 pixels and express the alignment profile of the 
minor axes of the tidal fields as a function of i? — R- Figure El plots (|u(i?) • u(i? )| 2 ) a t 
z = as open squares with Rq = 8.7/i _1 Mpc. As can be seen, as R decreases, the strength 
of the alignment gradually decreases. 



It was ILee fc Pen! (120011 ) who have for the first time shown that the auto-correlations 
between the principal axes of the linear tidal fields smoothed on two different scales can be 
approxi mated as a squa re- root of the density correlation coefficient, £(-Rq — R), defined as 



see also lLee et al.ll2009l ) 



C(Ro-R) = 








J P(k)W(kR)W(kR )d 3 k, 


4 - 


J P(k)W 2 (kR)d 3 k, 


4 = 


J P(k)W 2 (kR )d 3 k. 



(5) 
(6) 
(7) 



Adopting the model of ILee fc Pen! (120011 ) and extending it to the nonlinear regime, we 
approximate the spatial correlations of the minor principal axes of the nonlinear tidal fields 
smoothed on the two scales of Rq and R as 

(\u(R)-u(R )\ 2 ) = ^(Rq-R), (9) 
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where the power-law index, /3, is an adjustable parameter. In the original work of iLee fc Pen 



( 1200 ll ). the linear power spectrum is used for the evaluation of the density correlation coef- 
ficient, £(Rq — R), through Equations flSD-flSl)- In our work, however, we redefine £(Ro — R) 
in terms of the nonlinear density power spectrum instead of the linear power spectrum to 
account for the fact that the alignment profiles of cluster galaxies lie in the nonlinear regime. 
Fo r the calcula t ion of the nonlinear power spectrum, we adopt the analytic formula provided 



bv lSaito et all (120081 ). 



Adjusting the value of (3, we fit the numerical results to Equation Q via the standard 
^'-minimization. The solid and dashed curves in Figure [6] represent the best-fit models for 
the cases that the nonlinear and linear density power spectrum are used for the calculation 
of £, respectively. Figure M shows that the alignments between the minor principal axes 
of the tidal fields on different scales are indeed well approximated as a power-law of the 
density correlation coefficient especially in the large AR = Rq — R section and that the 
density correlation coefficient expressed in terms of the nonlinear power spectrum yields a 
better fit to the numerical results. The discrepancy between the numerical results and the 
best-fit analytic curves in the small- AR section (AR < 4 ft, -1 Mpc) should be attributed to 
the low-resolution of the Millennium tidal fields (256 3 grids on a periodic box of linear size 
500/i- x Mpc). 

As the alignments between the minor principal axes of the tidal fields on two different 
scales are now found to scale as a power-law of the density correlation coefficient, we model 
the three dimensional alignment profile of cluster galaxies as 

(caB 3 0)-± = Ae(AR), (10) 

where the amplitude A and the power-law index n are two adjustable parameters, quantifying 
how strongly the satellites are aligned with the external tidal fields at the virial radius and 
how fast the alignment profile drops with AR, respectively. Adjusting the values of A and 
n, we fit the numerical results presented in §2 to Equation ( ITOl) at each redshift. For the 
two dimensional case, the constant term in Equation (llOp is changed from 1/3 to 1/2. For 
the fitting model, we use the density correlation coefficient, £(AR), expressed in terms of 
the nonlinear density power spectrum and set the value of Rq at the maximum top-hat scale 
found from the clusters belonging to each sample. 

The best-fit values of A a nd n are determined through minimizing the following gener- 



alized x 2 (IHartlap et all 120071 1: 



x 2 



A cos2 0i C^j 1 ^ cos 2 6j, (11) 
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where A cos 2 Qi denotes the difference between the numerical result and the analytic model 
(Eq.JTU]) at the z-th bin, and (Cy) is the covariance matrix defined as = ((cos 2 ^ — 
(cos 2 9i)) (cos 2 9j — (cos 2 8j))), where the ensemble average is taken over the eight Jackknife 
resamples (see §2). The value of Rq for the best-fit analytic model is set at the maximum 
top-hat scale calculated from the masses of the clusters belonging to each sample. 

The best-fit amplitude and power-law index for the three and two dimensional cases 
are listed in Tables [TJ121 The uncertainties in the measurement of the best-fit parameter 
yalues represent the margina l ized errors computed as the curvature of x 2 a ^ its minimum 



(IBevington &: Robinsonlll996l ; I Wall fc Jenkins! 120081 ). Figures [2J15] plot the best-fit models as 



a function of 1 — R/Rq in the left panels, which reveal that our analytic models (Eq.[TD]) are 
indeed in good agreements with the numerical results for all cases. We also plot the same 
alignment profiles versus £(AR) in the right panels of Figures [2H to show explicitly that the 
alignment profiles scale as a power-law of the density correlation coefficient. 

The fitting results also show qualitatively how the alignment profiles of cluster galaxies 
change with mass scale and redshift. As the cluster mass scale increases, the amplitude of the 
alignment profile increases slightly while its power-law index decreases significantly. As the 
redshift increases, the amplitude of the alignment profile shows little change but its power- 
law index increases considerably. This result indicates that the satellite galaxies located in 
the massive clusters at lower redshifts tend to be more strongly aligned on average with the 
external tidal fields. To test the statistical significance of the dependence of the alignment 
profile on logM and z, we analyze a simple linear regression of the best-fit values of A and 
n listed in Tables [Q12] for the three subsamples: (I, II, III) and (I,IV,V), respectively. Table 
[3] which lists the slope, a, of each linear model shows quantitatively that the dependences 
of n on logM and z are both significant with \a\ > 1 while the dependences of A on logM 
and z are insignificant. 

The strong dependence of n on z can be explained as follows. At higher redshifts, 
the relative dominance of the dark matter density parameter, Q m (z) is greater and thus 
the nonlinear effect inside the clusters would turn on earlier, resulting in higher degree of 
the deviation of the major axes of the pseudo inertia tensors of the cluster galaxies from 
the minor principal axes of the external tidal fields. As for the strong dependence of n on 
log M, the more massive clusters are believed to form more recently and thus their relaxation 
process has yet to be completed. In consequence, the satellite galaxies of the more massive 
clusters would retain better memory of the preferential alignments with the minor principal 
axes of the external tidal fields. Meanwhile, the weak dependence of A on z and logM may 
be interpreted as follows. The amplitudes of the alignment profiles are largely determined 
by the degree of the alignments of those satellite galaxies located at the virial radii of their 
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host clusters. Since those satellite are likely to have just merged into the clusters, they 
should always exhibit the strongest alignments with the external tidal fields, regardless of 
the redshift and mass scale. 

It is worth mentioning here that our model for the alignment profile of cluster galaxies 
does not require precisely accurate measurement of the cluster's virial radii, given that during 
the fitting procedure the value of Rq for our model is set at its maximum top-hat scale from 
the clusters belonging to each sample. In other words, the alignment profile of cluster galaxies 
can be constructed even when we have just approximate measurements of virial masses of 
the sample clusters. 



4. A COSMOLOGICAL IMPLICATION 

Although we have considered only a single cosmological model to measure the alignment 
profile of cluster galaxies that lies in the non-linear regime, the strong redshift dependence of 
the alignment profile shown in §3 hints that it might depend on the background cosmology. 
Furthermore, the fact that the alignment profile scales as a power-law of the density correla- 
tion coefficient, £(AR), suggests that it might be useful for breaking the fi m -o"s degeneracy 
since £(Ai?) is independent of the power spectrum normalization o~ 8 . 

To explore this possibility, we first examine how the density correlation coefficient, 
£(AR), changes with the matter density parameter. Varying the value of Q m from 0.1 to 
0.5 and fixing the other key cosmological parameters at the values used for the Millennium 
simulations, we repeatedly calculate £ (AR) expressed in terms of the nonlinear power spec- 
trum (see §3), at z = 0, the results of which are plotted in the left panel of Figure [0 It 
reveals that the higher the value of Q m is, the more rapidly the correlation coefficient £(AR) 
drops with AR. This result is consistent with our finding in §3 that the alignment profile 
of cluster galaxies drops faster with AR at higher redshifts when the relative dominance of 
dark matter is high. 

The other key cosmologic al parameter that we pay attention to is the neutrino mass, m u . 



Recently, iMarulli et al.l (120111 ) have demonstrated by hydrodynamic simulations that if the 
neutrinos have non-zero mass (in a ACDM+z/ cosmology), then their free-streaming would 
have a strong impact on the formation and evolution of the large-scale structure, suppressing 
the formation of small- scale obj ects and thus decreasing th e abundance of galaxy clusters 



and their clustering strength (see iLesgourgues fc Pastorll2006l for a recent review). Yet, they 



emphasized that these effect caused by the neutrino free streaming are strongly degenerate 
with the power spectrum normalization factor, cr 8 . 



-li- 



lt is expected that if there are larger amount of free-streaming neutrinos, then they 
would defer the merging of the satellite galaxies into the clusters and thus the cluster galaxies 
would keep better the memory of the external tidal fields. The higher the fraction of massive 
neutrinos is, the higher amplitude and lower power-law index the alignment profile of cluster 
galaxies has. Under the speculation that our analytic model for a ACDM cosmology can be 
extended to A + z/CDM cosmology, we calculate £ (AR) for five different cases of the neutrino 
mass fraction, /„ = Q u /Q m . Regarding the nonlinear matter density power spe ctrum for a 
ACDM+z/ cosmology, we use the analytic formula provided by lSaito et al.l (120081 ). The right 
panel of Figure [7] plots the matter density correlation coefficient at z = for five different 
cases of /„. The higher the value of f u is, the correlation coefficient, £(AR), drops less 
rapidly with AR, consistent with our expectation. 



5. SUMMARY AND CONCLUSION 

We have determined the three and two dimensional alignment profiles of cluster galaxies 
using the halo and semi-analytic galaxy catalogs from the Millennium simulation. The 
alignment profiles of cluster galaxies are defined as the average alignments between the major 
axes of the pseudo inertia tensors of the cluster galaxies on two different top-hat scales Rq 
and R (where Rq is the top-hat radius enclosing cluster's virial mass and R < Rq). The 
alignment profiles of cluster galaxies reflect how well the satellite galaxies, after merging into 
the clusters, retain the memory of the preferential alignments with the minor principal axes 
of the external tidal fields and how fast they gradually lose the alignment tendency under 
the influence of the internal tidal fields as the cluster's relaxation proceeds and nonlinear 
effect becomes stronger. 

It has been shown that the alignment profiles of cluster galaxies expressed as a function 
of AR = Rq — R behave as a power-law of the density correlation coefficient. The alignment 
profiles exhibit weak dependence of the amplitude A on mass scale and redshift but strong 
dependence of the power-law index n whose statistical significance is tested by using a 
simple linear regression analysis. As the cluster mass decreases and the redshift increases, 
the alignment profile drops faster with AR, which is consistent with the scenario that the 
faster merging process of the satellites and earlier onset of the nonlinear effect inside the 
clusters tend to break more easily the preferential alignments of the cluster galaxies with 
the minor principal axes of the external tidal fields. The weak dependence of A on mass 
scale and redshift implies that the satellite galaxies located at the virial radii of their host 
clusters always keep good memory of the external tidal fields regardless of the mass scale 
and redshift. 



We have suggested that the alignment profiles of cluster galaxies might be a compli- 
mentary cosmological probe, being useful especially for breaking the fl m -<Js and f v -<T$ de- 
generacies (where f v is the neutrino mass fraction) owing to the insensitivity of £(AR) to 
cr 8 . To prove this theoretical concept, however, it will be necessary to test the validity of 
our model against N-body simulations for different cosmologies other than the ACDM model 
assumed for the Millennium simulations. As a large catalog of the galaxy cluster s and their 



satell ites from the Sloan Digital Sky Survey Data Release 7 is in the pipeline ( jHao et al. 



20101 ) . our future work is in the direction of conducting numerical and observational tests of 



the alignment profiles of cluster galaxies presented here. 
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Fig. 1. — Illustration of the anisotropic spatial distributions of the satellite galaxies in the 
two dimensional projected plane and the angle 9 between the major axes of the pseudo 
inertia tensors of the satellites within two different radial distances. 
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Fig. 2. — Three dimensional alignment profiles of satellite galaxies in cluster halos identified 
at z = from the three different mass bins: 14.0 < log(M / [h~ l M@\) < 14.2 (triangles); 
14.2 < log(M/[/i -1 M ]) < 14.4 (diamonds); 14.4 < \og(M /[h~ l M ]) < 14.6 (squares). The 
solid lines correspond to the analytic models, Equation f fTOi) . with the best-fit parameters. 
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Fig. 3. — Three dimensional alignment profiles of satellite galaxies in clusters with loga- 
rithmic mass in the range of 14.0 < \og(M/[h- l M & }) < 14.2 at three different redshifts: 
z = 0, 0.5, 1 (open squares, triangles and diamonds, respectively). The solid lines corre- 
spond to the analytic models, Equation (fTOj) . with the best-fit parameters at z — 0, 0.5 and 
1, respectively. 
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Fig. 4. — Same as Figure [2] but for the two dimensional case. 
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Fig. 5. — Same as Figure [3] but for the two dimensional case. 
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R -R (Mpc/h) 

Fig. 6. — Alignments between the minor principal axes of the Millennium tidal fields 
smoothed on two different top-hat scales (square dots) at z — 0. The analytic models 
expressed as a power law of the density correlation coefficients expressed in terms of linear 
and nonlinear density power spectrum (dashed and solid line, respectively). 
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Fig. 7. — The density correlation coefficient expressed in terms of the nonlinear density 
power spectrum for five different values of fl m and of f v (left and right panel, respectively). 
The other key cosmological parameters are fixed at the same values as for the Millennium 
simulations. 
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Table 1. Redshift, logarithmic mass range, number of the selected clusters, best-fit 
amplitude, and best-fit power-law index of the three dimensional alignment profiles of 
satellite galaxies in the clusters from the five samples 



sample 


z 


\og{M/[h- l M Q }) 


N c 


A 


n 


I 





[14.0, 14.2] 


1250 


0.570 ±0.006 


3.130 ±0.038 


II 





[14.2, 14.4] 


649 


0.600 ±0.010 


2.350 ±0.017 


III 





[14.4, 14.6] 


337 


0.630 ±0.006 


2.280 ±0.090 


IV 


0.5 


[14.0, 14.2] 


832 


0.580 ±0.005 


3.420 ± 0.041 


V 


1 


[14.0, 14.2] 


369 


0.560 ±0.005 


4.700 ± 0.042 
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Table 2. Same as Tabled] but for the two dimensional alignment profiles 



sample 


N c 


A 


n 


I 


1077 


0.470 ± 0.004 


2.930 ±0.049 


II 


552 


0.490 ±0.001 


2.190 ±0.021 


III 


289 


0.470 ± 0.004 


1.810 ±0.096 


IV 


724 


0.470 ± 0.003 


3.060 ±0.065 


V 


324 


0.480 ± 0.005 


3.840 ±0.050 
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Table 3. Result of the linear regression test of the dependence of the best-fit parameters 
of the alignment profiles on the logarithmic mass and redshift 



function 


a 


A(logM) 


+0.100 


ra(log M) 


-1.417 


A(z) 


-0.007 


n{z) 


+ 1.047 



